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^ '. Abstract 

Q ' We re-examine the model of Altarelli, Cabibbo, Corbo, Maiani and Mar- 

tinelli for inclusive semileptonic B decay, in the light of recent calculations in 
^ I heavy quark effective theory. The model can be shown to have no l/rrif, cor- 

rections, with a suitable definition of the b quark mass mi,. However, we find 
that the structure of the 1/m^ terms is incompatible with the predictions of 
heavy quark effective theory. The numerical significance of this discrepancy 
is discussed. 
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I. INTRODUCTION 



Recently, there has been much progress in the apphcation of heavy quark effective theory 
(HQET) to the calculation of inclusive semileptonic decays of B hadrons. It was shown by 



from deep inelastic scattering. In particular, one can perform an operator product expansion 
on the hadronic tensor contributing to Hi, — > Xiug, where Hi, represents a (heavy) bottom 
hadron, and X a sum over final hadron states. At leading order in this expansion, the parton 
model prediction results, while higher-order terms, which are suppressed by inverse powers 
of the b quark mass m^, give non-perturbative corrections. In addition, it was demonstrated 
the expansion is free of 1/ terms, essentially because in the effective theory, the operators 
with the correct dimensions to contribute to such terms have matrix elements which vanish 
to leading order. 

The 1/ml corrections to the parton model for the inclusive differential rate dT/dEi 
(where Ei is the final lepton energy in the B hadron center of mass) were calculated by 
Bigi et ai 0. The non-perturbative corrections enter the calculation in the form of matrix 
elements of higher dimension operators: the two that contribute in this case are 



where Hb{v) denotes a B hadron with velocity v, D is the QCD covariant derivative. Gap is 
the QCD field strength tensor and denotes the quark field in the heavy quark effective 
theory. 

Subsequently, a very detailed exposition of the 1 / corrections was given by Manohar 
and Wise |^. This paper presented double differential distributions (which also appeared 
around the same time in Q), and also considered polarisation effects in Af, decays. In 
addition, it was also shown that some of the higher order terms could be interpreted as 
resulting from averaging the zeroth-order (i.e. parton model) decay rate over the residual 
4-momentum of the bound b quark inside the B hadron. 



Chay, Georgi and Grinstein |l| that the problem can be attacked using methods familiar 




and 
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This interpretation is reminiscent of the older calculation of Altarelli et al 0, which 
we henceforth refer to as the the ACCMM model. In this model, bound state corrections 
are incorporated by introducing a wavefunction for the two quarks inside the B meson, and 
then averaging the partonic decay rate over their relative 3-momentum. The size of the 
ACCMM wavefunction is determined by the Fermi momentum pp of the quarks inside the 
B meson, which is of order Aqcd- Consequently, one can perform a large expansion of 
the ACCMM result, and investigate whether it can, with suitable choice of parameters, be 
made to agree with the QCD predictions of This procedure is the main focus of this 

letter. Throughout, we ignore the corrections of perturbative QCD to the process, discussing 
bound state corrections only. In Section 0, we introduce our notation and briefly review 
the averaging of Manohar and Wise. The ideas behind the ACCMM model are introduced 



in Section [II, where we also attempt to find a correspondence with the newer calculation 



to order We present our conclusions in Section |^ 



II. AVERAGING IN HQET 



From the charged- current interaction in the Standard Model 

V,b qjl^PLb h^^PLVt + h.c. 



V2 

■^v,,f;j,, + h.c. 



(1) 



(where j is a flavour index, V is the CKM quark mixing matrix, is the Fermi weak 
decay constant and Pl is the projection operator (1 — 7^)/2) one can derive the inclusive 
differential semileptonic decay rate 

Here one has summed over all final states Xj containing a quark of flavour j = u or c. The 4- 
vector q is just the sum of the 4-momenta p and p' of the final state lepton and antineutrino, 
while Ee and Ei, denote the energies of these particles. Throughout, we neglect the lepton 
mass rrii. The leptonic tensor Lajj is given by 



where our convention for the totally antisymmetric tensor is 60123 = +1- Wai3 can be 
expressed as 

W^p = E (Ph, -Px-q) (h, 4 (0) X,) (X, J^^ (0) H,) . (4) 



In this expression, we have normalised the B hadron state \Hb) as in ref. to remove 
extra factors of mass from the decay rate (2). In the parton model, Wa/s is approximated 
by its value from free quark decay: 

= 6 [ml- 2mb {E, + E,) + - ) (5) 
X { -Iga/B {mb -En- E^) + nibVaVp - ^ieaf^f^uV'^q"} ■ 

Here we have introduced the velocity of Hf,; in addition we explicitly note the dependence 
on rrib and because these will later be averaged over. As shown in reference 0], most of 
the I /ml corrections to the parton model result to Wa/s can be obtained by averaging the 
RHS of equation (5) over the residual 4-momentum of the bound b quark. If one defines 
the mass m[, velocity v'^ and momentum of the moving b quark via 

= m'bv'^ = rribVf, + ; (v'"^ = = ij , (6) 

then one can expand its contribution to Wa/s to second order in k: 



(7) 



Note that this expression contains an expansion of the 6 function in (5), which cannot be 
justified everywhere in the Dalitz plot. In fact, the OPE contains an analagous expansion 
(of the j-quark propagator), which breaks down in some of the regions where a partonic 
description is not reasonable 0. We will assume throughout this paper that we are in a 
region of phase space where this truncation is valid. 
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It turns out that almost all the terms coming from the OPE analysis up to order 1/m^ 
can be reproduced by replacing Wa/3 by the average (over k) of equation (7), with the average 
values 

(A;") = EbTTibv'^ (8a) 
= -^m^ (c,"^- vV) . (8b) 

Prom this point of view, the order of magnitude of the quantities Ei, and Kf, might not be 
obvious; however, they can be expressed as matrix elements of operators between B meson 
states. As a result, they are expected to be of order {Aqcb / iTiby . Incidentally, one should 
note that the corrections which cannot be accounted for via this averaging are proportional 
to the parameter Gh, which is defined in HQET by 



Gb = ^^{Hb{v,s) KgC^pa'^H, Hb{v,s)), (9) 



where s and v denote the spin and velocity of the heavy meson Hj,, and Zf, is a renormalisation 
constant. To leading order in l/rrih, this operator is a spin-spin interaction between the heavy 
quark and its surrounding "brown muck" . A simple physical interpretation of its effects has 
not yet been given, however. 



III. COMPARISON WITH ACCMM 

The ACCMM model introduces bound state corrections to the parton result for B meson 

decay by considering the meson to be made up of a b and spectator quarks. Purthermore, 
the energy of the spectator is assumed to be given in terms of its 3-momentum p^p by 

Esp = Vp'p + (10) 

where nisp is a free parameter of order Aqcd- The 4-momentum P'^ of the b quark is now 
fixed by the requirement that the sum of the b and spectator momenta be that of the 
B hadron. One finds that 
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Mb - E.. 



sp 



Mb- 



(11) 



P 



Psp 



and hence that the square of the invariant mass of the moving b quark is just 




M| + ml - 



2MB^P% + ml 



(12) 



To perform the averaging over the motion of the quarks inside the hadron, one intro- 
duces a wavefunction ipiPsp)- Typically, one expects that this function will have a width 
of order Aqcd and that it will fall off rapidly for values of momentum p^p larger than this 
value. Also, because we understand the B meson to be an L=0 state, we assume that V'(Psp) 
is spherically symmetric. In practice, one takes 



where iV is a normalisation factor and pp, the Fermi momentum of the quarks inside the 
hadron, is taken to be an adjustable parameter of order Aqcd- One can now write down the 
hadronic tensor in the model, by assuming it is just the average over the quark 3-momentum 
of the boosted free quark value: 



In this expression, we have explicitly included a Lorentz time-dilation factor in the integral, 
which was omitted in 0. It should also be noted that in the integral is cut off at the 
maximum value of \psp\ allowed in the decay; however, the wavefunction is so small at this 
point (the exponential factor in (13) is tiny there) that we can ignore the effect of this upper 
limit in what follows. 

Since the wavefunction factor in (14) effectively restricts \psp\ to be at most of order pp, 
we can envisage expanding the factor wj^^ {m'f^, v') / v'q in (14) by treating p^p, pp and m^p 
as small. Clearly, the first term in this expansion will just be W^^] {Mb,v). The fact that 
this depends on Mb and not ruf, is not accidental: the ACCMM model was specifically 
constructed to try to circumvent the dependence of the partonic decay rate on the unknown 




(13) 




(14) 
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parameter m^. Just as was the case in HQET, our expansion will contain average values, 
which are now given (as averages over p^p) by 

(/> = I d'pspmPsp)f fiPsp). (15) 

If the two types of averaging were to produce identical results to this order, we would 
be able to express the HQET parameters nib, Eb and Kb in terms of Mb, Pf and rusp- In 
fact, there would be some dependence on the ACCMM wavefunction through, for example, a 
quantity like (psp)- Because our expansion starts at Mb, and because the HQET expansion 
contains no l/mb terms, we will evidently have to absorb 1/Mb terms into the definition of 
the quark mass rrib. Evidently, in ACCMM we have from (11) 

{P") - V- {Mb - {^P% + m%)) ■ (16) 
Consistency with (6) and (8a) can only be obtained if 

rrib (1 + Eb) = Mb- (^^JT^) • (17) 
The 4-momentum k in (6) can now be expressed in the B rest frame as 

= [EbMB - ^Jp% + m% + (^^pl + ml'j , -p,,) , (18) 
so that we have, to leading nonvanishing order 




On the RHS of this equation, the first term comes from the 00-component in the B rest 
frame, while the second is the ij-th, taking into account rotational symmetry. The Oi- and 
iO-components can easily be seen to vanish. 

Since the two tensors in (19) are linearly independent, the only way to reconcile this 
equation with (8b) is to have 



and also for the coefficient of v^'v^ to be of order AQQj-,/m6, i.e. to vanish to the order we are 
considering. However, for a general wavefunction of size ppi this quantity will be of order 
In addition, one can compute the contribution of the extra term to the quantity LajsW""^ 
in (2), to make sure it does not vanish. We find a contribution to this quantity of 

2AmlE, {5\X) {QrribE, - AEj - AE,E, - q') 

+ 2 6"{X) {mt -Ee- E,f (2m,E, - q")] (21) 
with X = ml- 2mfc {Ei + E^) + - m^, 

where we have defined the dimensionless parameter by 

AbUil = (jplp + m%) + . (22) 

Incidentally, one should note that this contribution cannot just cancel the Lorentz time- 
dilation factor mentioned after (14), as that gives an overall factor (1 — Eb) to the decay. 
Consequently, it seems to be impossible to find an exact correspondence between the aver- 
aging of the phenomenological ACCMM model and the QCD calculations of HQET. 

IV. CONCLUSIONS 

We have analysed, in two different approaches, the corrections to the parton model for 
inclusive semileptonic B meson decay up to order l/ml. Although we have found that the 
ACCMM model cannot be put into exact correspondence with the HQET calculation, we 
should point out that the discrepancy between the two is likely to be small numerically. 
This is mainly due to the fact that l/m^ corrections vanish in both cases if equation (17) 
is satisfied. As demonstrated in 0, the I /ml terms adjust the parton result by around a 
percent in the region where the expansion makes sense. Also, the additional term appearing 
in (cf. equation (19)) in the ACCMM case is likely to be suppressed compared to 

the other term for a monotonically decreasing wavefunction. This term would dominate the 
other corrections in the m — oo limit, but for the B it is practically indistinguishable from 
the higher order terms we have ignored. 
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The ACCMM model's disagreement can be traced to its assumption that the spectator 
quark in the B meson can be treated as an on-shell particle of mass rrisp, of the order of 
the QCD scale. This led to one having to define rrib in such a way so as to avoid l/rrib 
terms in the expansion. However, making this definition was not sufficient to ensure that 
the 1/ml terms had the correct form. In the context of HQET, it is clear that the "spectator 
quark" of ACCMM corresponds to "brown muck" . However, it seems unlikely that the latter 
can be treated simply as an object of mass m^p. If one wanted to improve the model to 
make it compatible with HQET, one could start by relaxing the on-shell condition for the 
spectator quark. 

Finally, it is worth noting that equations (6) and (8), together with (18) and the way 
we defined the quark mass in ACCMM, suggest that, as far as these inclusive processes 
are concerned, the parameter Ef, can be absorbed into a redefinition of mf,. This can be 
confirmed by examining the expressions in section 6 of reference 0: the Eb correction terms 
have the form 

d 

Ebrrib— — (parton model quantity), 

onib 

where the partial derivative is taken at fixed lepton momenta. Presumably, this should still 
hold true in the case of r decay modes, where the lepton masses cannot be ignored . This 
observation might be useful if one were trying to fit the HQET parameters to an experimental 
spectrum, for example. 
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